(19) 



J) 



(12) 



Europaisches Patentamt |||||| 
European Patent Office 
Office europeen des brevets (1 1 ) EP 0 795 985 A2 

EUROPEAN PATENT APPLICATION 



(43) Date of publication: 

17.09.1997 Bulletin 1997/38 

(21) Application number: 97103953.1 

(22) Date of filing: 10.03.1997 


(51) IntCI. 6 : H04L 27/26 


(84) Designated Contracting States: 


(71) Applicant: NOKIA TECHNOLOGY GmbH 


DE FR GB IT 


75175 Pforzheim (DE) 


(30) Priority: 13.03.1996 Fl 961164 


(72) Inventor: Rinne, Jukka, 




c/o Nokia Technology GmbH 




75175 Pforzheim (DE) 



(54) Compensation of phase errors in multicarrier signals 



(57) The invention relates to a method for compen- 
sating channel errors in digital data communication in a 
signal in sampled data form. The method evaluates the 
phase error of the examined sample on the basis of the 



phase errors of at least two preceding samples, the esti- 
mate being used to compensate the phase error of the 
examined sample (y m (n)). 
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Description 

The invention relates to the method claimed in the preamble of claim 1 for compensating channel errors, such as 
errors, phase noise and frequency offset caused by-multipath propagation in a digital data communication system, 
especially an OFDM system (Orthogonal Frequency Division Multiplexing). 

In this disclosure channel stands for a transmission path used in data transmission, which usually is a radiopath in 
wireless data transmission systems. 

In wireless airborne data transmission, the signal will always be distorted to some extent in practice when passing 
from the transmitter to the receiver. Distortions are caused among other things by buildings, woods and meteorological 
variations. Moreover, as the receiver and/or the transmitter moves, the conditions will change constantly, thus affecting 
flawless signal reception. In practice, a signal reaches the receiver along a number of paths, and thus the distances 
covered by the signals following different paths will not be equally long, generating phase and amplitude differences 
between the signals. Imperfect oscillators also entail phase errors, for instance. 

The OFDM system uses several subcarrier wave frequencies, allowing several signals to be transmitted at the 
same time. Each signal x m serves to modulate one subcarrier wave. The modulated subcarrier waves are combined 
and transformed into a mode suitable to be transmitted on the transmission path. In the receiver, demodulation is car- 
ried out, separating the transmitted signals to various receiving channels. One known OFDM demodulator can be 
formed by using Inverse Discrete Fourier Transformation (IDFT), and accordingly, an OFDM demodulator can be 
formed with Discrete Fourier Transformation (DFT). Figure 1 illustrates the principle of the OFDM system as a sche- 
matic block diagram. In figure 1 the channel block represents a transmission path, such as a radiopath, and then h(n) 
represents n transfer functions of the channel. 

Prior art methods for compensating linear channel distortions are based on linear equalizers. In the OFDM system 
linear equalisation is most frequently performed in the frequency domain owing to the nature of the OFDM system. Prior 
art methods are mainly based on minimizing the mean square error (MSE) of received and sent samples. To achieve 
this, adaptive algorithms are used, such as the least mean square (LMS) principle. Another option is to use a zero forc- 
ing (ZF) criterion in the frequency domain combined with a proportional algorithm (PA), which in practice provides faster 
convergence and also easier selection of the learning constant than adopting the MSE criterion in connection with an 
LMS algorithm. 

Both the previously known methods mentioned above use the same type of method for compensating received and 
demultiplexed samples. The received samples are demultiplexed with Discrete Fourier Transformation (DFT), i.e. by 
transforming the received samples from time domain to the frequency domain, and subsequently the samples are mul- 
tiplied with complex numbers, which are adjusted to the algorithm applied. 

A known method that improves the compensation of constant phase errors uses minimum mean square error to 
estimate phase offset. This enables the frequency offset to be calculated and tuned to the oscillators. 

The problem of LMS and PA algorithms is that when a^small error is aimed at, the algorithms converge slowly. Then 
the algorithms are unable to follow fast variations in the channel characteristics. This causes problems especially 
regarding relatively long OFDM symbols, in which a phase change may be notable even during one single symbol dura- 
tion. 

The purpose of the present invention is to eliminate the above drawbacks and to achieve a method and a device 
for compensating errors caused among others by multipath propagation, phase noise and frequency offset. The method 
of the invention is characterised by the features defined in the characterising clause of accompanying claim 1. 

The inventive idea is that a linear channel distortion is compensated by first adopting the minimum MSE criterion 
and the LSM algorithm. The mean phase rotation can be calculated on the basis of channel distort! on- compensated 
signals and estimates corresponding to real values. After this the calculated mean phase rotation is used for the two 
preceding symbols, and then the phase error can be estimated with linear approximation in the symbol examined. The 
phase error being known, the phase noise of the symbol can be compensated in the time domain by inverse phase rota- 
tion. 

The present method achieves appreciable advantages over prior art methods. With the use of the method of the 
invention, phase noise can be compensated far more effectively than with currently known methods. In addition, the 
method of the invention allows compensation of larger frequency offsets than known methods. The method of the inven- 
tion can be implemented fully digitally, preferably in an application program of the receiver, and then no changes to the 
high frequency end of the receiver will be called for. 

The invention will be described in greater detail below, with reference to the accompanying drawings, of which 

figure 1 illustrates the OFDM system as a basic block diagram, 

figure 2 illustrates the principle of compensating a linear phase distortion as a coordinate presentation and 



figure 3 illustrates a schematic block diagram of the device for compensating phase distortion in accordance with 
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the invention. 

Currently known compensating methods are not adequate for compensating errors arising in time-variant commu- 
nication channels. In such channels frequency offset and phase noise will occur. A slowly changing phase error causes 
a common carrier phase error, which can be estimated. The use of the phase error estimate in the two preceding sam- 
ples y m ] ,y m _ 2 allows a linear curve to be fitted to these and to determine the phase error on the basis of this in sample 
y m to be received, and thus linear components in the phase noise can be compensated. 

The method of the invention is examined below from a mathematical point of view exemplified by an ODFM signal 
x m (n) , which is conducted over a channel, in which a phase error and multipath propagation are generated in signal 
x m {n). The oscillator noise is marked with the term z m {n), and it is assumed that the oscillator noise phase atone varies 
and the amplitude is constant, i.e. 

z m (n)=e''* m{n) (formula 1) 

If the OFDM symbol x m (n) is relatively long, i.e. the phase varies notably during one symbol duration, the phase 
error must be compensated. 

In a time-invariant multipath channel the channel response h m {n) can be given in the form 

h m {n)^a(X)6{n-X) = h{n),a(^C (formula 2) 

where 




1 ? when a — 0 
0, othenvise, 



r is the number of multipath propagated components and C denotes a set of complex numbers. 

Assuming that the noise generated after downconverting is neglectable, the received samples y m (n) are then given 

by 

yjn)=(x m (nyh(n) + d m (n))z m (n) (formula 3) 

where * denotes convolution. 

The compensation of received samples y m (n) is done in two separate steps. In the first step, the linear part of the 
phase error is compensated in time domain with the time compensation term t m n , giving the phase compensated sam- 
ple 

■/m( /7 ) SB /m(' 7 )Ufi (formula 4) 

In the second step, compensation is done in the frequency domain with the frequency coefficient C m (k), the fre- 
quency compensated sample being given by the formula 

X m {k)=Y m (k)C m (k) (formulas) 

where 

Y m (k)=DFT[y m (n))=l£y m (n)e n "* IN (formulae) 

n=0 

The output of the equalizer is obtained on the basis of samples X m {k), i.e. with the decision function P[X m {k)}. 

As the value of the frequency compensated sample X m (k) given by formula 5 may lie between the symbol values 
adopted in the system, the next corresponding value is selected among the symbol values available. This is known as 
such in discrete systems. 

The frequency compensating coefficients C m {k) are iterativefy solved with an LMS algorithm to meet the following 
criterion: 
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(formula 7) 



The updating rule then has the form 



(formula 8) 



W Tm nrv' 5 3 ' 6amin9 C °" Stant Wh ' Ch ' S 3 p0Sitive real number and ^variable during updatin 0 
15 z m W- p [X m W]-X m {k) and Y m *{k) is the complex conjugate of Y m (k). The learning constant A is selected such 
that the term in formula (8) converges. 

The linear phase error is compensated with the aid of the MMSE phase error compensation estimates of the phase 
errors e m0 and e m _ 2 of the two preceding symbols, i.e. 



20 



0 m k=0 



(formula 9) 
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Formula (9) enables the phase compensation term Q m to be solved, giving: 



e m =arg 



AM 



N-1 

Ei* m (*>i ! 

k=0 



(formula 10) 



where X m (k) is a demultiplexed, phase compensated sample equalised in the frequency domain and X m (k) is a 
demultiplexed, non-phase compensated sample equalised in the frequency domain. 

The required phase compensation for the nth sample of the mth symbol can be expressed as 



45 



where 



and 



f m [n)=a m n+b. 



/ m (-(3A//2 + 2/V g ))=e m . 2 



(formula 11) 



(formula 12) 



(formula 13) 



so where /V is the number of samples in guard interval. In other words, the phase correction is positioned at the centre of 
each OFDM symbol, and thus constants a and b can be solved from the above equations as follows: 



e m-r e m-2 



(formula 14) 



55 and 



,^3e m . r e m . 2 h2/v q (2e m , r e m , P ) 



2(N+N g ) 



(formula 15) 
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The estimates of the phase correction angles % m . 2 and-e m .-, have been defined such that offsets greater than n are 
converted into corresponding 2n complements. Figure 2 shows the principle of estimating phase offset. 

The coefficient required in the phase compensation of the received nth sample of the mth symbol can now be 
expressed as follows: 

f mn= e' Un, as e*"e tam '\ n=0...W-1 (formula 16) 

In practical implementation this is best expressed with the following formula: 

n 

t mn - e' bv Y{ e ' 3m (formula 17) 



is Here N complex multiplying operations and two e ,x operations will be necessary. 

Figure 3 shows a schematic block diagram of coefficient 1 of the invention. The received nth sample ym(n) of the 
mXh symbol is conducted through the IN line of correction coefficient 1 to the first Discrete Fourier Transformer 2 and to 
the linear phase compensator 3. The Fourier-transformed signal Y m (k) is conducted to the first linear equalizer 4, 
where the Fourier-transformed signal Y m {k) is multiplied with coefficient C m (k) as in formula (5) so that the term X m (k) 

20 is equalled by the term X m (k) and the term Y m (k).\s equalled by the term Y m (k). The output of the first linear equalizer 
4 provides sampled signals. X m (k).to be taken to the phase error calculator 5. 

In the linear phase compensator 3 the phase of an input signal Y m (n) is compensated by means of the correction 
coefficients a m , b m calculated in the phase error calculator 5 in accordance with the above formulas (4), (1 7). The signal 
y m (n) phase compensated in the time domain is conducted to a secondary Discrete Fourier Transformer 6, where the 

25 signal is transformed to the frequency domain. The signal Y m (k), transformed to the frequency domain and phase com- 
pensated, is taken to a second linear equalizer 7. In the second linear equalizer 7 the signal Y m {k) is multiplied with the 
coefficient C m (k) as in formula (5). 

The signal X m (n) provided by the second linear equalizer 7 is taken both to the decision member 8 and to the cor- 
rection coefficient updating member 9. The decision member 8 selects the system symbol corresponding next to the 

30 signal X m (n), this symbol thus corresponding to the output OUT of the equalizer 1. The selected symbol is also taken 
to the correction coefficient updating member 9 and the phase error calculator 5. 

The invention can also be implemented such that the secondary Discrete Fourier Transformation Y m {k) (formula 
6) is calculated for some subcarrier waves only. In this way, however, the phase error estimate will be less accurate and 
more subsceptible to noise, but in low noise cases the accuracy will be satisfactory. 

35 Instead of time domain, phase noise compensation can be executed in frequency domain using convolution. In this 
case, the secondary Discrete Fourier Transformation 6 is replaced with a digital filter, thus allowing a less complex 
structure for equalizer 1. 

The invention is not restricted to the above embodiments, but can be varied within the scope of the accompanying 
claims. 

40 

Claims 

1 . A method for compensating channel errors in digital data communication in a signal in sampled data form, charac- 
terised in that the phase error of the sample examined is estimated at least on the basis of the phase errors of the 

45 two preceding samples, the estimate being used to compensate the phase error of the sample (y m {n)) examined. 

2. Method as claimed in claim 1,' characterised in that the phase error estimation is performed in the frequency 
domain. 

so 3. Method as claimed in claim 2, characterised in that the phase error compensation is performed in the time domain 
with inverse phase rotation. 

4. Method as claimed in claim 1 , characterised in that the compensation of the phase error of the sample examined 
is performed in two steps, the linear part of the phase error being compensated in the time domain in the first step, 

55 and compensation being performed in the frequency domain in the second step. 

5, A method as claimed in claim 4, characterised in that in the first step compensation is done with the equation 
y m (n) =y m {n)t mn (formula 4) . where 
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Y m (n) = the nth phase compensated value of the mth symbol, 
y m {n) = the received nth sample of the mth symbol, 
trn n = the time equalising term, 

and in the second step compensation is performed with the equation X m (k)=Y m {k)C m (k) (formula 5), where 

X m (k) = the nth frequency compensated value of the mth symbol, 

Y m (k) = the Discrete Fourier Transformed sample of the phase compensated value (y m (n), 

C m {k) = the frequency coefficient. 

6. Method as claimed in claim 4 or 5, characterised in that the compensated value (P[X m (k)]) of the sample (y m (n)) 
examined is formed by selecting the symbol value available in the system corresponding next to the frequency com- 
pensated value (X m {k)) provided in the second step. 

7. Method as claimed in claim 6 P in which a guard interval (t g ) is added to the symbol (x m {n)) during transmission, 
characterised in that the time equalising term (t m n ) is calculated with the equation 



t m ^e if ^ n) =e ib »e ia ™ n ,n=0...N-1 



(formula 16) 



where 



' m-r u m-2 



(formula 14) 



^(3e m . r e m . 2 )+2A/ g (29 m , r e m , 2 ) 

2(N+N Q 



(formula 15) 



=arg 



/V-1 

ZP[*m(*)]*mW 

k*0 



Zi* m wr 



(formula 10) 



N g = the number of samples in guard interval (t g ), 
N - the number of samples of symbol (y m (n)). 

8. Method as claimed in claim 6 or 7, characterised in that the frequency coefficient {C m {k)} is iteratively solved by 
using the LMS algorithm to meet the following criterion: 



iV-i 



mm'Z\? m (k)CM-X m {kl( 



(formula 7), 



a frequency coefficient calculated with the following equation being used to compensate the subsequent sample in 
the frequency domain: 



" m+1 

(k)=CJk)+Ac m (k)Y m *(k) 



(formula 8) 



where 



A = learning constant lFT is invariable during updating and is selected such that the term in formula (8) con- 
verges, 

*mW = P[*m(A)]-* m (/t). and 
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Y m "{k) is the complex conjugate of Y m (k). 

9. Device (1) for compensating channel errors in digital data communication in a signal in sampled data form, char- 
acterised in that the device (1) comprises: 

5 

- means (2, 5, 6, 8) for estimating the phase error in the examined sample (y m (n) t y m {n)) on the basis of the 
phase errors of at least two preceding samples, and 

- means (3, 4, 7, 9) for compensating the phase error of the examined sample (y m (n), y m {n)). 

w 10. Device (1) as claimed in claim 9, characterised in that the means (2, 4, 5, 6, 8) for estimating the phase error of 
the examined sample (y m ( n ). ym( n )) comprise means (2) for transforming the examined sample {y m (n). y m (n)) to 
the frequency domain. 
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